土中爆炸成腔问题中确定炸点状态的一种计算反求方法

提出了一种通过给定的土中爆炸成腔毁伤效应确定炸点状态的计算反求方法。该方法将确定炸点状态的反问题转化为求解爆炸毁伤效应的计算值与给定值误差函数最小的优化问题。在反求过程中,采用基于误差减小比率技术的多项式近似模型代替土中爆炸数值分析模型,以便提高反求效率。采用Tikhonov正则化方法克服反求过程中出现的病态问题。在此基础上,引入信赖域管理策略判断当前近似模型与实际模型的逼近程度,以确定最优的反求向量。炸点状态反求结果与实验结果的对比分析表明,该方法能够有效且稳定地通过给定的毁伤效应实现炸点状态的反求,这可为炸点状态的设计提供参考。     

Variational assimilation in combination with the regularization method for sea level pressure retrieval from QuikSCAT scatterometer data I: Theoretical frame construction

A new method of constructing a sea level pressure field from satellite microwave scatterometer measurements is presented. It is based on variational assimilation in combination with a regularization method using geostrophic vorticity to construct a sea level pressure field from scatterometer data that are given in this paper, which offers a new idea for the application of scatterometer measurements. Firstly, the geostrophic vorticity from the scatterometer data is computed to construct the observation field, and the vorticity field in an area and the sea level pressure on the borders are assimilated. Secondly, the gradient of sea level pressure (semi-norm) is used as the stable functional to educe the adjoint system, the adjoint boundary condition and the gradient of the cost functional in which a weight parameter is introduced for the harmony of the system and the Tikhonov regularization techniques in inverse problem are used to overcome the ill-posedness of the assimilation. Finally, the iteration method of the sea level pressure field is developed.     

一种改进型多组分气体的Tikhonov正则化特征光谱提取方法

针对FTIR光谱分析数据计算量大,具有相同基团的多组分混合气体交叉灵敏度过高,以及在线分析过程中基线漂移等问题,提出了一种改进的TR2-1(Tikhonov 2norm-1norm)正则化特征变量提取法。该方法借鉴TR1和TR2模型的基本思想,引入谱线距离和谱线系数绝对值最小化两个约束项来保证所提取特征量的准确度,消除基线漂移所带来的影响,并结合上述两种模型建立了不适定问题的最优化求解通式。该通式有效地克服了经验法和实验法确定正则矩阵和参数所带来的误差。实验以浓度为0.01%~20%的甲烷、乙烷、丙烷、正丁烷、异丁烷、正戊烷和异戊烷气体为例进行了特征光谱选取。结果表明,对于浓度为0.2%的甲烷气体,其预测误差平方和仅为2.6×10-4,可决系数达到0.959 2,分析准确度高,有效增强了TR正则法的实用性。     

Multigrid Solution of a Lavrentiev-Regularized State-Constrained Parabolic Control Problem

A mesh-independent, robust, and accurate multigrid scheme to solve a linear state-constrained parabolic optimal control problem is presented. We first consider a Lavrentiev regularization of the state-constrained optimization problem. Then, a multigrid scheme is designed for the numerical solution of the regularized optimality system. Central to this scheme is the construction of an iterative pointwise smoother which can be formulated as a local semismooth Newton iteration. Results of numerical experiments and theoretical two-grid local Fourier analysis estimates demonstrate that the proposed scheme is able to solve parabolic state-constrained optimality systems with textbook multigrid efficiency.     

自适应非凸稀疏正则化下自适应光学系统加性噪声的去除

自适应光学系统可以实时测量并校正波前信息,但是系统中大量的噪声严重影响了系统的探测精度.自适应光学系统中一般为加性噪声,本文提出一种全新的变分处理模型去除加性噪声,该模型采用自适应非凸正则项.非凸正则项在保持图像细节上较凸正则项具有更好的效果,能更好地保持点源目标的完整性.另外,根据不同区域的噪声水平自适应地构建正则化参数,使不同区域的像素点受到不同程度的噪声抑制,可以更好地保持目标的边缘细节.在算法实现上,为了解决非凸正则项收敛性较差的缺陷,采用分裂Bregman算法及增广拉格朗日对偶算法进行计算.实验及数值仿真结果都表明,该方法能够较好地去除系统中的加性噪声,且光斑信号保存得较为完整,处理后的质心探测精度及信噪比较高.     

三维声场边界元法几乎奇异积分问题分析

以三维声场问题为例,提出一种准确计算高阶单元几乎奇异积分的半解析算法.首先分析高阶单元几何特征,构造近似几何量,然后应用扣除法,将奇异积分核函数分解为规则核函数与近似几何量表达的奇异核函数.规则核函数积分采用常规Gauss数值积分计算,奇异核函数积分采用半解析算法计算.给出三维声场内问题和外问题经典算例,计算了近边界点的声压,结果证明本文半解析算法的有效性和准确性.     

一种求解弹性l_2-l_q正则化问题的算法

给出了一种求解弹性l_{2}-l_{q}正则化问题的迭代重新加权l_{1}极小化算法, 并证明了由该算法产生的迭代序列是有界且渐进正则的. 对于任何有理数q\in(0,1), 基于一个代数的方法, 进一步证明了迭代重新加权l_{1}极小化算法收敛到弹性l_{2}-l_{q}(0相似文献   

Investigation of self-adaptive algebraic tomography for gas reconstruction in larger temperature range by multiple wavelengths absorption spectroscopy

We develop a self-adaptive algebraic tomography algorithm （SAATA） to investigate the simultaneons reconstruction of concentration and temperature distributions in larger temperature range from two views. The simplified optical arrangement with fewer projections is realized by extension of spectral information at multiple wavelengths, resulting in great potential in applications of practical combustion diagnosis. Tile results show SAATA can perform much better reconstructions in 300 3000 K temperature range than genetic simulated annealing algorithm and least-square orthogonal-triangular decomposition method with two- wavelength scheme. More phantoms are created to demonstrate the capability of SAATA to capture the peaks and adapt for both flat and sharp temperature distributions. Meanwhile, the advantage of high stability ensures better reconstruction performance at noise levels from 0.1% to 10% in projections.     

Casimir effect for a scalar field via Krein quantization

In this work, we present a rather simple method to study the Casimir effect on a spherical shell for a massless scalar field with Dirichlet boundary condition by applying the indefinite metric field (Krein) quantization technique. In this technique, the field operators are constructed from both negative and positive norm states. Having understood that negative norm states are un-physical, they are only used as a mathematical tool for renormalizing the theory and then one can get rid of them by imposing some proper physical conditions.     

初始模型对含噪动态光散射数据正则化反演结果的影响

分别采用最小模型矩阵、最平坦模型矩阵、最光滑模型矩阵作为初始化模型,对加入5种不同水平随机噪声的90nm窄单峰、90nm宽单峰和250nm窄单峰、250nm宽单峰颗粒体系的模拟分布进行了正则化反演,并对反演结果进行比较。结果表明:当噪声水平为0时,正则化初始模型的选择对反演结果没有明显影响。随着噪声水平的增加,采用三种初始化模型反演得到的峰值误差和粒度分布误差都随之变大,但采用最平坦模型和最光滑模型反演得到的峰值和粒度分布误差明显小于采用最小初始模型的反演误差。当噪声水平大于0.01时,选择最平坦初始模型获得的粒度分布结果优于采用最光滑初始模型和最小初始模型获得的结果,而采用最光滑初始模型反演得到的峰值优于最平坦初始模型和最小初始模型的反演峰值。因此,采用正则化算法处理含噪动态光散射数据时,为得到最优的粒度分布信息,宜采用最平坦初始模型,若需要获取最准确的峰值信息,则应选择最光滑初始模型。